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Abstract 

The absence of self averaging in mesoscopic systems is a consequence of long-range intensity 
correlation. Microwave measurements suggest and diagrammatic calculations confirm that the 
correlation function of the normalized intensity with displacement of the source and detector, 
AR and Ar, respectively, can be expressed as the sum of three terms, with distinctive spatial 
dependences. Each term involves only the sum or the product of the square of the field correlation 
function, F = Fg. The leading-order term is the product, the next term is proportional to the 
sum. The third term is proportional to [F(AR)F(Ar) + [F(AR) + F(Ar)] + 1]. 

PACS numbers: 41.20.Jb, 05.40.-a, 71.55.Jv 
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Short-range correlation in waves transmitted through random media is manifest in the 
intensity speckle pattern. The leading contribution, C\, to the cumulant correlation function 
C of intensity normalized to its ensemble average on the output surface of the sample is 
given by the square of the field correlation function, C\ = F|. Neglecting internal reflection 
from the surface, its dependence upon displacement Ar on the output surface is given by 
Ci(Ar) = F|(Ar) = F(Ar) = (sin kAr/kAr) 2 exp(— Ar/£ s ), where k is the wave vector and 
i s is the scattering mean free path |l| . This term dominates intensity fluctuations. Defining 
a correlation length, 5r, as the first zero of C\ gives 5r = n/k = A/2. The intensity is 
correlated far beyond 5r as a result of scattering within the medium || [| f§, ||, ||, [7| so that 
intensity values in remote speckle spots are not statistically independent. This gives rise to 
two additional contributions to C, which can therefore be expressed as, C = C\ + C 2 + C 3 
||, ||, and leads to greatly enhanced mesoscopic fluctuations 0. The C*2 term produces 
an enhancement in total transmission fluctuations over that given by the field factorization 
approximation by a factor of L/£, where I is the transport mean free path and L is the 
sample length || [7j. The C3 term is the source of universal conductance fluctuations, which 



are enhanced by a factor of (L/£) || 10 1. The magnitude of C\ at a point is unity as a result 



of normalization, whereas the magnitudes of Ci and C3 are expansions in 1/g with leading 
terms of order 1/g and 1/g 2 [[|, respectively, where g is the dimensionless conductance. 
Since the onset of localization is at g — 1 [|n|] , the two terms beyond the field factorization 
approximation for C, reflect the approach to localization. 

In this Letter, we use microwave measurements and diagrammatic calculations to show 
that each of the contributions to C may be expressed in terms of the square of the field 
correlation function with regard to displacements of the source, AR and detector, Ar. 
The C\ term is F(AR)F(Ar), the C2 term is proportional to [F(AR) + F(Ar)], while 
the C 3 term is proportional to [F(AR)F(Ar) + [F(AR) + F(Ar)} + 1]. When the in- 
tensity correlation is considered at a shifted frequency, Au, the full correlation function 
remains a sum of three terms each being a product of the corresponding terms in C, thus, 
Ci = Aj(Az/)Cj(Ar, AR) (i = 1,2,3). Absorption alters the magnitudes of C2 and C3, but 
it does not change the spatial structure of these terms. 

Initial measurements of angular intensity correlation, carried out in the far field of weakly 
scattering media, gave C, which was essentially equal to C\ |T2|, [TJ|. Recently, measurements 
of the spatial correlation of the field on the sample surface have yielded the C\ contribution 
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directly [|H |. Measurements of intensity correlation between points on the sample surface 



and the interior of the sample on a scale greater than the wavelength have allowed the 
observation of C2 || • In addition measurements have been made of the spectral correlation 
of the C\ [|rjj and C2 terms fEB] , |T7| , |T^] and of the temporal correlation of the C\ |ll|, C2 [2C 



and C3 PTJ terms in colloidal samples. However, the variation of C on a subwavelength scale, 



as well as the structure of correlation with displacement of both the source and detector, 
have not been reported previously. The present measurements and calculations allow us 
to discern the structure of intensity correlation and to relate it to the correlation of the 
underlying field. 

Measurements are made in a disordered dielectric sample contained within a reflecting 
tube. Radiation of frequency v emitted by a source antenna at R at one end of the tube is 
detected at a point r at the other end. We denote the intensity at r due to a source at R 
by I u (f,R) and consider the normalized cumulant correlation function, 

C(Ar, AR) = (5Ur-R)5I„ir',R'))/(I u (r,R))(I„ir\R')) , (1) 

where 5I V is the deviation of the intensity from its ensemble average value, Ar = \r — r'\ and 
AR = \R — R'\ are the displacements across the output and input surfaces, respectively, and 
(. . .) denotes the average over an ensemble of random realizations. The leading contribution 



to C obtained by factorizing the fields is M, 0, 0, 122 



d(Ar, AR) = [(^(^^^(f',^))! 2 /^^^))^^',^)) • (2) 

The samples studied are random mixtures of 1.27-cm-diameter polystyrene spheres at a 
volume filling fraction of 0.52. They are contained inside a 100 cm long copper tube with a 
diameter of 7.5 cm capped with thin Plexiglas face plates. Field and intensity spectra are 
obtained in the frequency range 16.8 - 17.8 GHz with a step size of 625 kHz using a vector 
network analyzer. Measurements are taken in an ensemble of 690 random samples by rotating 
the tube to create new configurations of the spheres after each set of spectra are taken. Field 
spectra are obtained in each sample realization by translating an antenna detector along a 
line to each of 50 locations separated by 1.06 mm on the output surface for each of two fixed 
antenna sources at the incident surface, which are separated by AR = 3 ±0.1 cm = d. Field 
spectra are also taken by translating the source along a line for each of two fixed antenna 
detectors at the output surface, separated by Ar = d. At the separation d, the short-range 



term has largely decayed. The antennas are aligned perpendicular to the line of translation. 
Intensity spectra are obtained by squaring the field spectra. 




10 20 30 40 50 



At (mm) 

FIG. 1: Plots of C(Ar, AR = 0) and C 1 (Ar,AR=0), and theoretical fit to C. 

The spatial variations of C(Ar, 0) and C\(Ar, 0) are shown in Fig. [TJ. The C\ contribution 
is directly obtained by squaring the field correlation function, shown in the inset of Fig. |T|. 
Subtracting C\ from the full intensity correlation function gives the difference [C — Ci], 
shown in Fig. |2] for a single source ( AR = 0) and for two sources separated by AR = d. This 
difference includes the terms beyond the field factorization approximation. Measurements 
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FIG. 2: plots of [C - C x ](Ar) for AR = and 3 cm, and Ci(Ar) at AR = 0. 

of C(0, AR) and Ci(0, AR), i.e. for a fixed detector and a scanning source, have also been 
performed. Within experimental error, C(Ar, 0) and C(0, AR) were found to be identical 
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functions of their respective arguments. The same is true for C\. Similarly, we find that 
plots of [C — Ci] versus AR with Ar = and Ar = d are nearly the same as those shown 
in Fig. 2. Thus AR and Ar can be interchanged as required by reciprocity. 
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FIG. 3: Comparison of Ci(Ar), normalized by its value at Ar = 0, for AR = 3 cm and Ci(Ar) at 
AR = 0. 

Measurements of C\(Ar, AR) for AR = and AR = d are presented in Fig. ||]. Within 
the noise level of 10~ 4 , the two functions have the same variation with Ar, C\(Ar,d) = 
2 x 10~ 3 Ci(Ar, 0). This numerical factor, Ci(0, d), is equal to the value of C\(d, 0) within 
the uncertainty in d. This result, taken together with the aforementioned symmetry with 
respect to interchanging Ar and AR, suggests that C\ can be written as the product of two 
identical functions, Ci(Ar, AR) = F(AR)F(Ar). 

We now examine [C — C*i](Ar, AR), which is do minated by C*2 in our sample. This 
function is seen in Fig. ^ to fall to nearly one half its value when either AR or Ar increase 
beyond Sr when there is no displacement of the other variable. This shows that C*2 is 
given by the addition of two equal terms. The comparison of the short-range variation of 
[C — C%] with G\ in Fig. || suggests that the additive form factors are identical to F, so 
that the dominant contribution to C — C\ is proportional to [F(AR) + F(Ar)]. This would 
imply that C2(Ar) approach a constant value for Ar > 5r, whereas the measurement of 
[C — Ci}(Ar,d) is seen in Fig. ^| to fall slightly with increasing displacement. This could 
be the consequence of a slight departure from a quasi-lD geometry at the output face of 
the sample. There, the average intensity is slightly larger at the center than at the edges 
since the wave can spread beyond the cross section of the tube. Notwithstanding this effect, 
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our experimental results suggest that both G\ and C2 can be expressed in terms of a single 
form factor F(x), where x stands for either Ar or AR. C\ and C2 contain, respectively, the 
product and the sum of two form factors. In Fig. 0, the correlation function [C — Ci](Ar, d) 
is seen to lie above the dotted curve, which is proportional to C\. This suggests a constant 
contribution to C . For Ar > 30 mm, the correlation function becomes negative, but here 
the noise becomes larger than the signal because of the reduced number of pairs of points 
with increasing Ar. Such a long-range correlation for large values of AR may be part of C3. 
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FIG. 4: Semilog representation of the spatial and frequency dependence of C\ (a) and [C — C±] (b) 
for AR = 0. 



The structure of the joint spatial and frequency dependences of C\ and C%, is obtained 
from measurements of the correlation functions Ci(Au, Ar) and [C — C\](Au, Ar) for AR = 
0, shown in Figs. |]a and [|b, respectively. The semilog representations in Fig. |] show that, 
within the limits set by the noise level, Cj have the same frequency dependence for any 
Ar, while Cj have the same spatial dependence for any Av for i = 1,2. Thus their spatial 
and spectral variations for a single source are given by Ci(Av, Ar) = Ai(Au)Ci(Ar). The 
noise level found in C\ is low compared to that in C2 because the field correlation function, 
Fe, is computed and then squared to obtain G\, giving a signal to noise ratio which is 
the square of that for the field correlation function. The form of the intensity correlation 
function suggested by experiment is borne out in diagrammatic calculations, which are 
briefly summarized below. 

The three terms in C may be represented diagrammatically. The diagram corresponding 
to the C\ term describes two non-interacting diffusons attached to pairs, GG*, of averaged 
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Green's functions |7|, |25| . This diagram factorizes into a product of two field functions: 



(E„(r,R)E„,(r',R')) = / d 3 ri d 3 r 2 G„(r, f^G^r', r x )T vv ,(r x , r 2 )G„(r 2 , R)G* u ,(r 2 , R') , (3) 



where T vv i{f\, f 2 ) denotes the diffusion ladder and integration is performed over f±, f 2 inside 
the tube. For the quasi one-dimensional geometry, the diffuson is independent of its trans- 
verse coordinates, whereas the Green's functions decay rapidly on a scale of the mean free 
path. Taking R = R' and v = z/, we obtain 

(E u (f,R)E u (f',R)} = (^fj J d^G^rJGlir'^JiW^R)} ee F E (Ar) , (4) 

where I u (r, R) is the intensity at f, normalized to its average value at the output face of the 
tube. Thus, C 1 (Ar, AR) = F%(Ar)F%(AR) = F(AR)F(Ar). 

The diagrams corresponding to the C 2 and C 3 terms describe two incoming and two 
outgoing diffusons which interact in the bulk of the medium. In these diagrams, each pair 
GG* of external Green's functions contributes a spatial form factor F& as in Eq. (|4j). These 
give 

C = C l + C 2 + C 3 

= Ai(Av, a)F(AR)F(Ar) + —A 2 (Av, a) \F(AR) + F(Ar)] 

+ -|_A 3 (Az/, a) [1 + F(AR) + F(Ar) + F(AR)F(Ar)] , (5) 

where g is the average conductance of the sample. The coefficients Ai {i = 1, 2, 3) depend on 
the absorption coefficient a and the frequency shift Au. The structure of Eq. (|^) is similar 
to that of the correlation in transmission, obtained in the multichannel formalism 0, [|, 0. 
In the present case, however, all terms are described by a single spatial form factor. The 
field factorization term, A\(Av = 0, a) is unity and independent of absorption by definition, 
while A 2 (Au = 0, a) is given by [|, ^3|, ||J]: 



3 

A 2 (Au = 0,a) 



16a 



sinh 2a — 2a (2 — cosh 2a) 



sinh 2 a 



(6) 



The coefficient A 3 (a) depends weakly on a and its limiting values are A 3 (0) = l,A 3 (oo) 
15/16. 



For our samples, g ~ 7 and a ~ 3 Ul4| . From the measurement of [C — C\\ at Ar = 0, 
-r-A 2 = 0.076 with A 2 = 0.87. This is in agreement with Eq. (@) and calculated corrections 
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due to localization effects || [R| |3|, |4], p5| . Using the measured Ci(Ar) as the functional 
form F(Ar), following Eq. (|5|) and neglecting C3, we obtain a good fit of the spatial structure 
of the intensity correlation function, as shown in Fig. [1|. 

These considerations have applications to present efforts to enhance the capacity of wire- 
less communication by utilizing multiple antennas to detect the multiply scattered field 
ptj| , Antenna separation should be larger than Sr and the number of statistically inde- 
pendent antennas equals the inverse of the degree of long-range intensity correlation, 3g/2. 

In conclusion, we have uncovered the connection between the field and intensity correla- 
tion functions in the spatial structure of the three contributions to C. In contrast to the case 
of angular correlation ||, intensity correlation can be expressed in terms of a single form 
factor obtained from the field correlation function. We have demonstrated the multiplicative 
character of C\ and the additive character of C2. Calculations predict a mixed character for 
C 3 , which includes a multiplicative, an additive, and a constant term of equal amplitude. We 
observe the infinite-range component of C3 in the residual correlation when both the source 
and detector are displaced by more than the correlation length. Determining the proper 
breakup of C into its components is of particular importance when considering simultane- 
ous variations in space, time and frequency. Each term is a product of the corresponding 
C\, C2 and C3 correlation function in the appropriate variables. 
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